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ABSTRACT 

We apply a likelihood analysis to the data of Lauer & Postman 1994 With P(h) 
Q\ parametrized by (<r 8 , T), the likelihood function peaks at <r 8 ~ 0.3, T <J 0.025, indicating 
at face value very strong large-scale power, though at a level incompatible with COBE. 
There is, however, a ridge of likelihood such that more conventional power spectra do not 
^ seem strongly disfavored. The likelihood calculated using as data only the components of 
the bulk flow solution peaks at higher <r 8 , in agreement with other analyses, but is rather 
broad. The likelihood incorporating both bulk flow and shear gives a different picture. 
^ The components of the shear are all low, and this pulls the peak to lower amplitudes as 
T^j- a compromise. 

The Lauer & Postman velocity data alone are therefore consistent with models with 
^ very strong large scale power which generates a large bulk flow, but the small shear 
(which also probes fairly large scales) requires that the power would have to be at 
very large scales, which is strongly disfavored by COBE. The velocity data also seem 
^ compatible with more conventional P(k) with 0.2 <^ T <^ 0.5, and the likelihood is peaked 
i around <r 8 ~ 1, in which case the bulk flow is a moderate, but not extreme, statistical 

^ fluctuation. 

^ Applying the same techniques to the data of Riess, Press, & Kirshner 1995, the 

results are quite different. The flow is not inconsistent with the microwave dipole and 
we derive only an upper limit to the amplitude of the power spectrum: <r 8 <J 1.5 at 
roughly 99%. 

Subject headings: Cosmology: theory - observation - galaxies: clustering - distances 
and redshifts 
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1. Introduction 

There have recently been several analyses of the 
large-scale velocity data for 119 Abell clusters within 
15000 km/s obtained by Lauer and Postman (1994, 
hereafter LP) which suggest that most currently pop- 
ular models of structure formation (e.g., CDM and its 
variants obtained by adding tilt, a cosmological con- 
stant, or an admixture of hot dark matter) could not 
produce the large magnitude of the bulk flow seen, at 
the "2-3cr level" (Tegmark, Bunn & Hu 1993, Strauss 
et al. 1994, Feldman & Watkins 1994). All of these 
have focused on the bulk flow statistic, a highly re- 
duced representation of the data. Here, we will per- 
form a likelihood analysis which uses more of the data. 
We take as our "hypothesis space" the family of power 
spectra P(k; ag, T), parametrized by the usual am- 
plitude and shape parameters ag and T (Efstathiou, 
Bond & White 1992); see Eq. (7)), and calculate the 
likelihood £(cr s , T) = _P(data|<78r). We also calculate 
the likelihood for the case when the data are reduced 
to either the bulk flow alone (as in other analyses) or 
the bulk flow and a gradient. Since the gradient or 
shear also probes large scales, we would expect these 
to give results compatible with the bulk flow alone. 

More recently, Riess, Press & Kirshner (1995, here- 
after RPK) have analyzed data from a nascent survey 
employing the light curves of Supernovae la as dis- 
tance indicators; although the sample of 13 galaxies 
out to 7000 km/s is quite sparse, the errors (« 5%) 
are small enough to make the data interesting. We 
apply our techniques to this dataset as well. 

2. The Likelihood Function 

The main statistical tool we will use to analyze the 
peculiar velocity data will be Bayes' Theorem, and 
through it, the likelihood function, as in Kaiser 1988. 
Bayes' theorem can be written as 



p(6\DI) 



P (6\I) P (D\6I) 
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where p(a\bc) roughly means "the probability [den- 
sity] of a given b and c." Here, 6 represents the param- 
eters of the theory we are considering (here, ag and 
T), D represents the data, and I represents any prior 
information we bring to the problem. Thus, p(0\I) is 
the infamous "prior distribution" of the parameters, 
p(D\6I) = C{9; D) is the likelihood function, and the 
remaining factor in the denominator merely serves as 
a normalizing constant so that J d6 p(6\DI) = 1. Eq. 



(1) tells us how we update the probabilities one as- 
signs to various hypothetical power spectra. Other 
information (e.g., the normalization from COBE — 
although we do not use this information in this paper) 
can go into the prior. 

The data are given by a list of particle positions, 
velocities, and errors, calculated from the data. LP 
use the brightest galaxy in a cluster as a distance 
indicator, with a phenomenological relationship be- 
tween the absolute magnitude, L m , and the slope of 
the brightness profile of the galaxy, a. RPK use the 
shape of the SNIa light curve as a predictor of its 
luminosity. 

We shall assume that the initial density and veloc- 
ity fields had a Gaussian distribution with some power 
spectrum P(k) = (|i5fc| 2 }. Assuming linear evolution, 
the velocity power spectrum is 



(2) 



for a universe with critical density; for SI ^ 1, this 
would be modified by an additional factor of /(SI) 2 ~ 
SI 1 ' 2 . Under the assumption of Gaussian initial con- 
ditions, these power spectra provide a complete sta- 
tistical description of the density and velocity fields. 
We assume that the errors are Gaussian, as well. We 
can therefore write the likelihood function as 
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where S m are the line-of-sight peculiar velocities for 
the m = 1 to N data points, and R mn = (S m S n ) 
defines the correlation matrix for the data. We can 
split this into two independent terms: a "theoretical" 
covariance matrix 

R(v) 

and an "error" matrix 
If the galaxies are located at positions r m = r m r m , 
the theoretical covariance matrix is given by 



AirPdk 
(2tt) 3 



Pv (k)fmn (&) • 



(4) 



where we sum over the spatial indices i, j, but not 
over the cluster indices m, n, and the function f mn (k) 
is the angle-average given by 
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The other part of the correlation matrix simply adds 
in the velocity errors: 

R ml = ( a m + vl)bmn no SUIT1 (6) 
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where <r m is the velocity error for cluster m. The 
term cr* is added to take into account a variety of ef- 
fects: LP report a 1-d dispersion of 271 km/s between 
the BCG redshift and the average cluster redshift; we 
also must take into account the possibility that the 
velocity field has undergone some nolinear evolution. 
For field galaxies, this effect may be large, but for 
clusters it is expected to be small. In total, we choose 
(7„ = 350 km/s for the results we present; however the 
results are not strongly dependent on this value — the 
approximately A6cz errors in the BCG velocity are 
dominant; for the SNIa, even the 5% distance indica- 
tor errors still dominate, although less so. We note in 
passing that we have not taken into account the er- 
rors in the positions of the clusters (which are equal 
to those in the velocities, since cz = H^r + r • v); 
for the LP data, we do, however, use the expres- 
sions from Colless (1995) to calculate the "Hubble 
redshift" czh ~ cz — S that the galaxy would have 
in the absence of peculiar motions, which is in turn 
used to calculate the comoving distance to the clus- 
ter d ?s czh/Ho- We are also using Colless' fit to the 
L m — a relation, which assumes a random distribution 
of velocities aside from the bulk flow. 

We model the spectrum as an initial (n = 1 
Harrison-Zel'dovich) power law times an appropriate 
transfer function: P(k) oc kT 2 (k) We choose the fol- 
lowing ansatz, after Efstathiou, Bond & White 1992: 

T(k) = [i + (ak/r + (bk/rf 2 + (cfc/r) 2 ) "] ~ 1,v 

(7) 

with a = 6 Ah' 1 Mpc, b = 3.0/1" 1 Mpc, c = U/i" 1 Mpc, 
v = 1.13, leaving two free parameters, the overall am- 
plitude of the spectrum, given for example by <7s, and 
the value of T. The latter controls the location of the 
turnover from the large-scale power law P(k) oc k to 
the small-scale P(k) oc k~ 3 . For a universe with cold 
dark matter, T ?a Q,h, and T can be also be related 
to quantities in universes with decaying neutrinos or 
mixed dark matter. 

We have checked the method with realizations of 
Gaussian density and velocity fields, related by linear 
gravitational evolution. We do not simulate a par- 
ticular distance indicator relation, but add unbiased 
Gaussian errors on the velocity realization, of compa- 
rable magnitude to those in the LP sample. We are 
able to recover the value of ag to better than "one 
sigma" (i.e., 68%), but T is less tightly constrained 
(as it is for the actual data below). 

In Fig. 1, we show the results of applying the 



method to the LP data. In order to help remove 
any possible nonlinear signal in the data, we have 
grouped the data using a friends-of-friends proce- 
dure with a linking scale of 2500 km/s. We per- 
form an unweighted average of the clumps, assum- 
ing uncorrelated gaussian errors, leaving about 50 
clumps of galaxies. Recall that throughout we use 
SI = 1; in general the results apply to the combination 
/(£!)<78 ~ Q 06 /bg, the usual (3 parameter. The maxi- 
mum likelihood value for the parameters is ag ~ 0.3, 
T ~ 0.025; the data actually prefer low amplitudes for 
us with a peak in the power spectrum at very large 
scales. At this point, we have \ 2 = SmSni?"^ = 47.5 
for 53 degrees of freedom, the number of clumped 
clusters. 

We show contours of constant likelihood, with like- 
lihood ratios relative to maximum of 0.250, 0.201, 
0.140, 0.079, and 0.019. Assuming a constant prior for 
p(ag\I) and p(T\I) over the plotted area, these likeli- 
hood ratios are such that the contours enclose 50%, 
68%, 90%, 95%, and 99% of the accumulated poste- 
rior probability. (In addition, we show the 68% con- 
fidence limits for the COBE quadrupole Q rms -ps = 
17lg fiK from Gorski et al. 1994, which can also be 
expressed as a normalization for the power spectrum 
amplitude, calculated on a much larger scale than as, 
or the effective window function W 2 (k).) 

Preference for a low T might have been anticipated 
from other analyses of the bulk flow, which seem to 
require a large amount of large-scale power. There 
is, however, a ridge in likelihood along increasing T. 
Canonical CDM, with ag = 1,T = 0.5, has a likeli- 
hood ratio of 0.311, placing it within the 50% contour. 
A model with ag = 1,T = 0.2 has a likelihood ratio 
of 0.376, also within the 50% contour. Note that we 
have assumed Gaussian errors; if the actual distribu- 
tion has more power in the tails, then any limits may 
be considerably weaker than those quoted here. 

[If we do not perform the clumping procedure, the 
results are even more restrictive, compressing the con- 
tours toward the ag = axis at all values of T. At 
T ~ 0.2, the unclumped data require <jg <J 1.6 at the 
99% contour. In essence, this is an upper limit on 
detected power; however, we expect these results to 
be dominated by close associations of clusters where 
our linear results may not be applicable (Croft & Ef- 
stathiou 1994).] 

Because LP (or Colless 1995) use the data itself to 
calibrate their BCG distance indicator, assuming that 
the flow is modelled by only a bulk flow, we must in- 
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vestigate possible biases induced in using the individ- 
ual peculiar velocities derived as residuals from BCG 
L m — a relation. We have done this by using the 
calibration derived from both the best fit bulk flow 
and that assuming the CMB rest frame; the results 
are insensitive to this change — the mean L m — a re- 
lation does not change substantially no matter what 
flow model is chosen. This applies to both the shapes 
of the likelihood contours and the values for the bulk 
flow and shear derived below. 

Now, we turn to the RPK data; with only thirteen 
points within 7000 km/s, there is not much leverage 
on the shape of the power spectrum (i.e., T), so we 
have simply set T = 0.5 and concentrate on a%. In 
Figure 2, we see that the RPK data in fact give an 
upper limit on ag; there is no clear detection of power, 
only consistency with moderate amplitudes for the 
power spectrum. 

It is perhaps surprising that the RPK data, with 
only 13 galaxies, can put such a strong constraint 
on the amplitude of the power spectrum; this can be 
traced to the small error (ps 5%) in the SNIa distance 
indicator, a factor of three better than that of LP, 
somewhat making up for the factor of ten in sample 
size. 

3. Bulk Flow and Shear Likelihoods 

Why is the evidence against CDM-like models so 
different from conclusions drawn from analyses of the 
bulk flow in the LP data? To answer this, we now per- 
form a likelihood analysis with reduced descriptions 
of the data, incorporating bulk flow and shear. In ad- 
dition to allowing a comparison with other results, we 
expect the bulk flow and shear to be less affected by 
nonlinear evolution than the full cluster/galaxy ve- 
locity field (Croft k Efstathiou 1994, Bahcall, Cen k 
Grammann 1994). 

We model the peculiar velocity field as «;(r) = 
Ui + rjpij + • • • where Ui is the bulk flow, pij the shear 
tensor, and the series could obviously be extended to 
higher moments as desired. Alternatively, we could 
model the full velocity field including the Hubble ex- 
pansion by the same bulk flow and an anisotropic ex- 
pansion term, + rjHij + • • • with a "Hubble tensor" 
Hij = Hobij +pij ■ However, we retain the trace of pij , 
because our estimators are different than those used 
by LP, and so we would not predict quite the same 
Hubble constant (although the difference is small; see 
Table 1). 



Of course, we only have line-of-sight velocities, so 
we really only model S(r) = r;V;(r) and we are only 
sensitive to the symmetric part of the shear tensor, 
as expected. Nonlinearity and observational error can 
be modeled as an additional error term, which we will 
take to be a gaussian. 

The maximum likelihood values for Ui and pij (and 
any higher moments desired) can be considered as 
members of a 9-component vector a a ; we are linearly 
fitting the S(r n ) to the nine independent functions r;, 
riVjr, which we shall write as g p (r), p = 1 . . .9, with 
coefficients a p . The maximum likelihood values for 
the a p are given by the usual linear fit for gaussian 
errors a p : 

r iv^ S n g q (r„) _ \ - g P (r n )g q (r n ) 

n n n n 

(Lynden-Bell et al. 1988, Kaiser 1991). This expres- 
sion can be converted to the integral 

a p = J d 3 r W ip (r)vi(r); 

n ^ ' 

(Restricted to p = 1 . . . 3 so g p = r p , this is equivalent 
to the expression in Kaiser 1988) Because the tensor 
window function is a linear filter on the field Vi(r), 
the vector a p obeys the same Gaussian distribution, 
with a suitably altered correlation matrix, 

/Ajrk 2 dk 
■^r^.W^WH"! 1 ( 10 ) 

with an angle-averaged squared window function given 

w 2 Pq ( k ) = / ^ m P {^)w jq ^)kh (ii) 

First performing this calculation for the three com- 
ponent (bulk flow only) fit, we show the results in 
Fig. 3, along with components of the bulk flow in 
Table 1, which are consistent with those found by 
LP. This confirms the results of other investigators 
who have found that the large magnitude of the bulk 
flow, pointed in a direction away from the principle 
axes of the error matrix Aij , supports a large am- 
plitude for the power spectrum, and places theories 
like CDM, or MDM, normalized to COBE, at approx- 
imately "three-sigma" — between the 95% and 99% 
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probability contours. Note that there is a ridge of 
high likelihood that continues to very low T and ag 
(i.e., with approximately constant CMB quadrupole 
Q). The bulk flow only tests the power spectrum on 
very large scales, causing this degeneracy between the 
parameters. 

In Figure 4, we show the likelihood for the LP data 
based on a nine-point fit to both bulk flow and shear, 
as well as the components of the fit in Table 1. This 
fit spreads the peak of the likelihood function down 
to lower amplitudes — it is not the very smallest scales 
which are forcing the high-amplitude favored by the 
bulk flow, but the intermediate scales probed by the 
shear. (This is analagous to the high "Cosmic Mach 
Number" of Ostriker and Suto (1990), who compare 
the magnitude of large-scale flows to the velocity dis- 
persion on small scales.) As discussed in Feldman 
& Watkins (1994), the bulk flow window function 
is peaked at scales k -1 ;> 100/i -1 Mpc, but has a 
non-negligeable tail extending to small scales. The 
shear window function (more precisely, the appropri- 
ate components of the combined bulk flow & shear 
window function) is peaked at 30/i -1 Mpc <J k -1 <J 
lOO/i" 1 Mpc, also with a significant tail extending to 
smaller scales. 

The maximum likelihood value of the shear has 
components of the same order as the diagonal terms 
in the error matrix (A^) 1 ! 2 ; the data are unable to 
strongly distinguish from the shear expected for no 
clustering power, P(k) = 0. Thus, including this in 
the likelihood calculation pulls down the implied am- 
plitude of the power spectrum — the detection of the 
absence of shear is not the absence of a detection of 
shear. Again, there is a ridge of high probability ex- 
tending to very low T and ag indicative of the still- 
large scales probed by the shear. 

Taking into account both the shear and bulk flow 
still favor a somewhat larger amplitude for the power 
spectrum at moderate values of T — even the small 
shear is not enough to compensate for the large bulk 
flow. Only when we consider all of the data, as above, 
is the amplitude consistent with <78 ~ 1; this implies 
that the LP flow is rather "cold" aside from the bulk 
flow. 

In Figure 2, we show the likelihood for the RPK 
data for the bulk flow and shear. In this case, the 
bulk flow alone shows a mild detection of power (since 
there is a bulk flow!), consistent to about one sigma 
with the CMB dipole, but the shear + bulk flow to- 



gether only give an upper limit, since the shear is 
small. Note that, because the sample is so sparse, the 
window functions are quite broad and have significant 
contributions on all scales, so the "shear" and "bulk 
flow" measured by RPK are quite different quanti- 
ties than those measured by LP (Feldman & Watkins 
1995). 

4. Discussion 

Several other analyses of the bulk flow calculated 
from the LP data set have been presented. LP them- 
selves, along with Strauss et al. 1994, conclude that 
its large magnitude is a significant problem within the 
framework of all currently-favored theories of struc- 
ture formation, when normalized to the COBE ob- 
servation; they rule out these conventional theories 
(many of which correspond to points on our ag-T 
plane) at 94-98% confidence limits using "frequentist" 
Monte Carlo techniques. Feldman & Watkins 1994 
use some of the same techniques as Sec. 3. above, sim- 
ilarly ruling out these theories at the 95% or greater 
level. Our results are similar; when the bulk flow 
alone is considered, conventional theories are indeed 
disfavored at 95-99%; once more information is taken 
into account, however, these theories become signifi- 
cantly less unlikely. 

These analyses have used "frequentist" statistics: 
they calculate some quantity for the LP data, such 
as the bulk flow or a \ 2 that takes into account the 
spatial distribution of the data. Then, they calculate 
the probability that a particular model (e.g., partic- 
ular values of ag and T) would produce a value of 
this statistic as large or larger than observed, given 
the known properties of the error distribution of the 
data. In some sense, these analyses integrate the like- 
lihood over possible data sets. We perform a Bayesian 
analysis which uses the same probability distribution 
(our likelihood function) but integrates only over pos- 
sible parameters as in Eq. 1. The two methods are 
not equivalent except in the case of linear models 
with uncorrelated errors, which this most definitely is 
not — the model, encoded in the quantity Rij, appears 
only through the correlation matrix. In any case, 
when analyzing the same data (i.e., the bulk flow), 
we stress that the two methods do agree — popular 
power spectra such as CDM are disfavored. Note that 
it is not immediately clear how to extend the frequen- 
tist analysis to the case with a small shear and large 
bulk flow — do we want the probability of observing a 
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smaller or larger value of some suitable x 2 ? 

Methodological differences aside, we agree with 
these analyses in the following sense: on the largest 
scales probed by the LP data, the perturbation am- 
plitude appears to be quite large, as implied by the 
large bulk flow. On the intermediate scales probed 
by the shear, however, the amplitude is somewhat 
smaller — most of the components of the shear ten- 
sor are within "one sigma" of a purely poisson veloc- 
ity field (i.e., they are of comparable magnitude to 
the diagonal components of the error matrix of 
Sec. 3.). In order to fit best on all scales at once, 
we are forced to quite reasonable values of the ampli- 
tude of the power spectrum. Note, however, that this 
does not address the possibility that this particular 
parametrization of the power spectrum is inadequate 
to describe the data well. 

For comparison, we observe a bulk flow U = 842 ± 
398 km/s and an rms shear p = (pijPij) = 12 ± 10 h 
km/s/Mpc, where the large errors come from the 
sparse sampling of LP's data (i.e., the large mag- 
nitude of the error part of the correlation matrix, 
R ( mne)). Standard CDM (a 8 = 1, T = 0.5), sam- 
pled with LP's window function, predicts U ~ 568 
km/s and p ~ 14km/s/Mpc. A power spectrum with 
(78 = 2, T = 0.5 that reproduces the observed bulk 
flow predicts a somewhat larger and non-negligeable 
p ~ 19km/s/Mpc. (None of these results have been 
corrected for LP's "error bias," and the errors all 
assume a diagonal correlation matrix between the 
points; for this reason we do not quote similar num- 
bers for the RPK sample where the off-diagonal cor- 
relations are non-negligeable.) 

Thus far, we have assumed that the errors on the 
line-of-sight velocities are simply diagonal; there are 
no correlations between them. This neglects the fact 
that the L m — a relation is calibrated from the data it- 
self, resulting in correlations between the line-of-sight 
velocities, since the standard candle luminosities de- 
pend on the whole data set. Moreover, the relation 
is calibrated assuming that the velocity field is well- 
modelled by a bulk flow superposed on essentially un- 
corrected velocities (or uncorellated luminosity flu- 
cutations, depending on the details of the analysis as 
in Colless (1995)). This could potentially bias the re- 
construction of any velocities other than the bulk flow 
(conversely, the presence of strong motions on smaller 
scales would bias the bulk flow result). As mentioned 
above, we have allowed the L m — a relation to vary, 
assuming different bulk flow models, and the likeli- 



hood contours do not change substantially. We also 
note that the size of the scatter in the L m — a relation 
is not an issue, even though it may in fact dominate 
over the actual velocities, as long as the distribution 
of the errors is well-understood. 

Ideally, we should work with the observed L m and 
a observations directly and perform our own calibra- 
tion, marginalizing over unknown nuisance parameters — 
like the paramaters of the L m — a relation itself, which 
are not well-determined. (We pause here to note that 
"nuisance parameter" is a technical term, referring 
to parameters involving physics in which we are not 
presently interested.) 

For the RPK data, the situation is less severe; they 
calibrate their distance indicator with a separate sam- 
ple of galaxies, so these biases should not occur, al- 
though the parameters of the distance indicator rela- 
tionship should in principle be marginalized over in 
this case as well. In any case, this is a statement 
about the technical ease of adapting the procedure 
of RPK to this context; ideally, all of the available 
information should be used in both cases. 

Finally, how, if at all, do we reconcile the LP and 
RPK observations? With a sample of 13 galaxies, 
RPK observe a flow which they claim has nearly "con- 
verged" to the CMB dipole at the 7000 km/s depth 
of the survey, whereas LP find a strong bulk flow ap- 
proximately 90° away from the CMB dipole. The two 
surveys seem to be probing quite different velocity 
fields, not surprising due to their very different ef- 
fective window functions. Obviously, the discrepancy 
between this result and that of LP needs to be better 
understood: if it is not merely a rare statistical fluc- 
tuation, we may trace it back to some un-accounted 
for systematic bias in one of the distance indicators. 

The authors would like to thank Dick Bond and 
Michael Strauss for valuable discussion and especially 
thank Marc Postman and Tod Lauer for providing 
their data as well as insight into its use. We also wish 
to thank Matthew Colless for providing the fits to his 
BCG distance indicator relation. 
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Fig. 1. — Full Likelihood as a function of power spec- 
trum parameters <jg, and T for the clumped LP data. 
Contours of constant likelihood are plotted (see text 
for likelihood ratios relative to maximum). For uni- 
form priors over the plotted area, the contours enclose 
50%, 68%, 90%, 95% and 99% of the accumulated 
probability. Also shown are 68% confidence intervals 
for the COBE quadrupole (solid lines originating at 
the origin). We have assumed £1 = 1; otherwise the 
horizontal axis should be replaced by the combination 
/(fi)<7 8 = Cl° 6 a a = /3. 



This 2-column preprint was prepared with the AAS L^T^X 
macros v3.0. 
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Fig. 2. — Likelihood function for RPK data. The 
three curves correspond to the full likelihood, the bulk 
flow + shear and the bulk flow alone. Because they 
use progressively less information, the full likelihood 
is the most restrictive, followed by the bulk flow + 
shear and the bulk flow alone. 
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Fig. 3. — Bulk Flow Likelihood for the unclumped LP 
data, reduced to the 3 points of the bulk flow velocity. 
Contours of constant likelihood are shown, enclosing 
the same integrated probability as in Fig. 1. 
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Fig. 4. — Bulk Flow + Shear Likelihood for the un- 
clumped LP data, reduced to the 9 points of the bulk 
velocity and gradient. Contours of constant likelihood 
are shown, enclosing the same integrated probability 
as in Fig. 1. 
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Bulk Flow [km/s] 




Shear [h km/s/Mpc] 




U X Uy U Z 


Pxx 


Pxy Pyy Pxz Pyz 


Pzz 



Bulk 

± 

Bulk & Shear 
± 



530 
301 
484 
312 



80.0 
327 

-216 
364 



649 
232 
644 
241 



-6.96 
6.10 



-5.18 
4.36 



-1.48 
7.17 



-2.49 
2.88 



2.76 
3.48 



-1.22 
3.71 



Table 1: Components of bulk flow (relative to the CMB) and shear for the LP dataset. The rows labeled "±" refer 



to the diagonal components of the error matrix \I A P q 
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